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Abstract 

This paper proposes an agent-based model which reproduces different 
structures of animal groups. The shape and structure of the group is the 
effect of simple interaction rules among individuals: each animal deploys 
itself depending on the position of a limited number of close group mates. 
The proposed model is shown to produce clustered formations, as well as 
lines and V-like formations. The key factors which trigger the onset of 
different patterns are argued to be the relative strength of attraction and 
repulsion forces and, most important, the anisotropy in their application. 
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1 Introduction 

Group behavior in animals has greatly interested scientists and researchers in 
the past, and has received further attention in the last decades as a test case of 
self-organization. Recently it has attracted attention not only from biologists 
and ethologists, but also from physicists, mathematicians, and engineers. This 
interest has produced a huge amount of literature, which is well documented 
and reviewed. See, for instance, [521IM1II3] and [1]. Loosely speaking, the basic 
idea behind these works is that complex collective behavior arises from simple 
interactions among close animals. Following this idea, the aim of this paper 
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is to investigate a simple model of interactions within groups which is able to 
reproduce rather different patterns and structures. 

Let us briefly introduce the main features of our model and its relationships 
with literature. A detailed description is postponed to Section [5] First, our 
model is agent-based, in the sense that each animal is singularly considered. 
Furthermore it is leaderless, meaning that all animals act following the same 
set of rules and their behavior is not imposed by others. These assumptions 
are widely considered in literature, and accepted as biologically suitable for a 
variety of species. See, among others, ^Tl [5T1 [51 [7] . 

Second, our model is purely based on attraction-repulsion interactions be- 
tween group mates. Attraction allows the group to be formed and stay tight, 
while repulsion allows to avoid collisions between group mates and keeps them 
well spaced. In the majority of papers, attraction and repulsion are combined 
with velocity alignment. Here we keep aside the issue of alignment which has 
received a considerable attention in itself [40l [9] , and focus on (the superposition 
of) attraction and repulsion. Our approach is close to [HI in this respect. 

Third, each animal interacts with a limited number of group mates. The idea 
of having a limited number of interacting neighbors is not new. The work [18] 
already considers attraction towards the nearest neighbor, while later experi- 
mental investigations found interaction with the closest two- four individuals [T] . 
This fact has been included in several models, among others [ITl [5T1 . 

but it has not always been included in recent models [1^1 [HI [Ml H] , in favor of 
a purely metric notion of neighborhood: interactions occur among group mates 
which are less than a threshold apart. Very recently, the former idea has again 
been brought to attention by [2] , where the authors present experimental results 
regarding fairly large flocks of starlings. Results show that interaction occurs 
with up to six-seven neighbors, no matter how far they are. The same paper 
also gives some simulation results, which suggest that a better cohesion of the 
flock can be guaranteed in that way. We also want to point out a third ap- 
proach to neighborhood deflnition, based on Voronoi partitions. Its application 
in biology dates back to [TH], and it is well documented in the physicists [TS] 
and engineers [3] literatures. 

Fourth, each animal only interacts with group mates which are in a suitably 
defined sensitivity zone. Restricting the interactions to a sensitivity zone raises 
the issue of defining its size and shape. On this matter, there is a significant 
amount of works considering limited visual (or sensing) fields. A limitation in 
the animal's angle of vision has been incorporated in most models, assuming 
a blind rear zone [HI [331 HI ttH HH] : oi distinguishing between front and rear 
sensitivity [16]. In [24] authors assume repulsion and alignment regions to be 
elliptical (taking into account body shapes), and include blind areas. In the 
recent paper |30| , the authors clearly distinguish between visual field and sensi- 
tivity zones, stating that the behavioral rules they use in the model apply only 
in the front zone. However, up to our knowledge, the anisotropy of sensitivity 
zones has always been taken just as a given constraint, and not as a potential 
resource able to shape the group geometry. Here is the main contribution of 
our paper: showing that a restricted sensitivity angle can be a key element 
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in determining the structure and shape of an animal group. Indeed, we show 
that changing two parameters, namely two sensitivity angles for attraction and 
repulsion, we obtain cluster formations, line formations and V-like formations. 
Note that similar patterns have already been obtained by means of other mathe- 
matical models (see for example [THl HSl [30] ) , the novelty here is that our model 
is able to reproduce all of them, depending on few parameters. We believe that 
this can offer new biological insights. 

Besides investigating the main issue of the role of anisotropy, we also present 
some results on the effects of limiting the number of considered neighbors, and 
on the influence of the relative strength of attraction and repulsion on the inter- 
animal distance. The latter problem relates to the results in 

The observations drawn from simulations are accompanied by a formal math- 
ematical investigation of the model. One can find that a popular approach to 
the analysis of agent-based flocking models consists in defining a suitable poten- 
tial function, called "virtual" or "artificial" potential, whose gradient gives the 
dynamics. This variational approach has been followed by both engineers [38j 
and biologists [IHl US leading to important results. However, it is not possible 
to apply it to our model. With respect to other models, ours has two main 
features: state- dependent switching dynamics, and asymmetry of interactions. 
The former is due to the state-dependent definition of the set of group mates 
which interact with a given animal. This discontinuous state dependence has 
been included in previous literature: see for instance the treatment in |38| . 
based on non-smooth potentials. The latter feature is due to the fact that 
the limited number of interacting neighbors and the restricted sensitivity an- 
gles imply that interactions need not to be symmetrical (reciprocal). This fact 
prevents the application of the virtual potential approach. The paper [34] is a 
partial attempt in this direction, because it considers a second-order system in 
which velocity alignment is achieved by asymmetric interactions. Nevertheless 
agents' positions are controlled using symmetrical information exchange. As a 
consequence, its results are not useful to us, as we are interested in spatial con- 
figurations. Overall, we conclude that the mathematical analysis of our model 
requires a novel approach, which has to include discontinuous and asymmetric 
interactions among animals. 

The rest of the paper is organized as follows. The details about the model are 
given in Section[2j Then, extensive simulation results are presented in Section[3l 
whereas SectionHcontains the mathematical analysis. Later, Section[5]discusses 
the implications of our findings and their biological soundness. We conclude 
presenting some lines of future research. 

2 Model definition 

The animals in the model are represented by point particles, which have simple 
continuous-time dynamics. Given N €N, for alH € {1, . . . , N} and t G M>0; let 
Xi (t) € represent the position of the i-th animal, whose evolution is described 
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by the differential equation 



i,it) =v,ix{t)), (1) 

where x{t) is the vector {xi{t), . . . ,XN{t))- As in [23, we use a coordinate 
system moving with the group centroid. This means that we are modeUing 
relative movements of the individuals and group structure, rather than its global 
motion. The velocity Vi{x) is the sum of two contributions, expressing the effects 
of attraction and repulsion, 

v,ix) = vt{x)-vl{x). 

In more detail, each of these contributions depends on the relative position of 
the other animals, 

vlix) = ^ M\\x, 

The following definitions have been used. 

• The function /a : M>o — > M>o (resp., /r : R>o — > IR>o) describes how each 
animal is attracted (resp., repelled) by a neighbor at a given distance, 
assuming [| • || denotes the Euclidean norm in R^. 

• The attraction neighborhood (resp., the repulsion neighborhood 7?.") 
is the set of the n animals closest to the «-th one, which are inside the 
attraction (resp., repulsion) sensitivity zone. 

The above model is very general, and we need to specialize it by choosing 
the interaction functions and the shape of the sensitivity zones. We make the 
following assumptions. 

Al. The functions /a and /j. are assumed to be 

M\\Xj -X^\\) = F^\\Xj -X.,\\ , fr{\\Xj - X^\\) = -, 

\\Xj Xt\\ 

where -Fa and are two positive constants. 

A2. The sensitivity zones are depicted in Figured] and illustrated as follows. 
Let the center point be the animal's position, and let the horizontal axis 
(arrow-headed) represent the direction of motion. Attraction is active in 
a frontal cone whose width is given by the angle eta S (0, 360°] (dashed 
line). Repulsion is active both inside a disk of radius Rsr > {short-range 
repulsion) and in a frontal cone of width € (0°, 360°] (solid line). We 
stress that the sensitivity zones do not necessarily coincide with the visual 
field of the animal. They rather represent the zones which attraction and 
repulsion are focused on. 



1 1 \ -^i Xi 
\\Xj -Xi\\ 



\\Xj -x^W 
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Figure 1: The shape of the sensitivity zones. 



A3. The speed of each animal \\vi\\ is bounded from above by a constant v, 
The above assumptions resuh in the system 

Xi{t) = ^li ^ {Xj 

Some remarks are in order. 

Rl. The definition of the interaction neighborhoods Af and TZf aUows to have 
a priori bound on the number of effective neighbors, and therefore on the 
sensing and "computational" effort which is required for each animal. This 
fact, which copes with animals' intrinsic limitations, has been experimen- 
tally observed in biology, for fish [T] and birds [5]. The latter paper calls 
this neighborhood definition topological, as opposed to metric definitions, 
based on distance only. 

R2. Our assumption of unbounded sensitivity regions docs not intend to imply 
that animals sensing capabilities extend on an unlimited range, but rather 
that group dynamics happen in a relatively small area. 

R3. If Q!a = ttr = 360°, the parameter i?sr has no effect, and the interaction is 
completely isotropic as in the simulations presented in [5J[S]. Note that, 
even if there is no preference for any specific direction, the limitation of 
the number of considered neighbors makes the interactions not reciprocal, 
i.e. the fact that the i-th animal interacts with the j-th does not imply 
that the j-th interacts with the i-th. 



(2) 
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R4. By specializing the functions /a and /r, one can obtain various interac- 
tion models. Indeed many proposals can be found in the literature, as 
reviewed in [JT] and [H]. However, the shape of these functions is not 
the main point in our paper, and thus we have decided to focus on a sim- 
ple choice, in order to highlight the innovative part of our approach, i.e. 
the angle- dependent interactions. Similar considerations are valid also for 
some features introduced in other models, such as a neutral zone around 
animals [321 [S] or a hierarchical decision tree which allows the repulsion 
force to have the priority over the attraction force [HI [S] . 

3 Simulations results 

In this section we make use of the agent-based model ([2]) in order to show 
the effects of the anisotropic interactions on the shape of the group. Namely, 
depending on the angles a^, aj., and the ratio between repulsive and cohesive 
forces, we shall obtain either clusters, or lines or V-like formations. These 
patterns are described in the sequel. 

To perform the simulations, it is instrumental the introduction of the con- 
stant ^ = which allows to rewrite ^ as 

In equation ([3]), the unit of length is chosen to be the body length (BL) of 
the animal, and the time unit (TU) is the inverse of i^a- Simulations arc ob- 
tained solving the system of equations ([3]) via an explicit forward adaptive Euler 
scheme. Each run starts from a randomly generated initial configuration (con- 
tained in square of edge L), and ends when the system reaches a steady state. 
To take into account the uncertainties in sensing and motion of the animals, 
we include small additive random disturbances on the direction of the velocity, 
uniformly distributed in [— anoiso, cknoisc]- AH the steady-state configurations de- 
scribed in the sequel are robust to such noise. A summary of the parameters 
and their values is given in Tabic [TJ In the sequel we discuss the role of A^, n, 
^, tta and a^, which are most significant to us, whereas i?sr, Umax, Q^noisc and L 
are kept fixed. 

3.1 Clusters 

In this paragraph, we describe simulation results when interactions are assumed 
to be isotropic. These results are not dissimilar from others in the literature 
(e.g., [351 US])- ■we include them for two reasons. First, for comparison with 
the less usual patterns described in the following paragraphs. Second, because 
they allow some interesting remarks about the role of the model's parameters ^ 
and n. 
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Table 1: Model parameters. 



Name 


Symbol 


Unit 


Values explored 


Forces ratio 


? 


BL 


0.1-20 


Attraction angle 




degrees (°) 


0-360 


Repulsion angle 




degrees (°) 


0-360 


Number of animals 


N 


adimensional 


2-200 


Number of considered neighbors 


n 


adimensional 


1-iN-l) 


Short-range repulsion radius 




BL 


1 


Maximum speed 


^max 


BL/TU 


2-30 


Size of initial domain 


L 


BL 


15 


Noise magnitude 


^noisc 


degrees (°) 


0-10 



Thus, let us assume that a^^ = = 360°. As a consequence, the outcome of 
the simulations is a cohesive and well spaced cluster. For a better understanding, 
we make use of an indicator which is largely used in the literature (see for 
instance [211 [24]): the mean distance to the nearest neighbor ISilSiD , defined as 

1 ^ 

NND = -^min||x.-rE,||. 

We investigate the dependence of NND on the parameter ^: simulation results 
are shown in Figure [2l We observe that NND is an increasing function of ^, in 



N=30 N=100 




Figure 2: NND as a function of the ratio ^, for different values of n. Error bars 
denote variance across individuals. Plots assume = 30 (left) and A^ = 100 
(right). 

particular it increases roughly linearly in ^. This linear dependence is observed 
for any choice of n. Two other features are noticeable: first, if n = 1, animals 
asymptotically converge to a comfortable distance which is equal to ^. Second, 
for any fixed ^, NND decreases as n increases. Moreover, all these remarks do 
not depend on N. 
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These results can be compared with those in [25]: in that paper, the au- 
thors assume that all animals in the group interact among each other, and they 
conclude that, the larger the group, the closer packed it is. Our simulations, 
instead, suggest that the significant parameter is n, the number of neighbors 
which is taken into account, rather than N, the global number of animals. How- 
ever, from the biological point of view, the number of neighbors n is not truly a 
free parameter: n can not be too large because of the limited sensing and anal- 
ysis capabilities of the animals, and can not be too small either. For instance, 
a low value of n sounds unsafe from the point of view of collision avoidance. 

Figured] also shows that the variance of the distance to the nearest neighbor 
among the animals is quite small. This means that groups are internally uni- 
form, in terms of spacing. This uniformity is also apparent if one looks at the 
steady-state configurations: three examples are shown in Figure[3j Nevertheless, 



(a) Crystal-like cluster: n = 1 (b) Disordered cluster: n = 7 (c) Round cluster: n = N — 1 

Figure 3: Examples of internal structures of clusters obtained with ^ = 10 and 
N = 60. The groups arc moving horizontally from left to right. 

as n varies, significant modifications are apparent in terms of relative positions. 
If n = 1, the animals are deployed to form an hexagonal lattice, reminiscent of 
a crystal. For intermediate ri's the internal structure is rather disordered, and if 
n = TV — 1 it is made of concentric circles. Note that crystal-like patterns have 
also been found in |15j and in |25j , where they are compared with the less regu- 
lar structures obtained in [6] . Such patterns might be of interest for engineering 
application to environmental deployment of robots [4] or sensors [12] , 

3.2 Lines 

In this paragraph we show how restricting the sensitivity field (i.e. reducing 
and aa) induces the formation of an elongated group. The key element for the 
formation of these patterns is a restricted frontal sensitivity field. Let us denote 
by e the oriented elongation of the group (sec for instance |161 124|). defined as 
the ratio of the vertical to the horizontal side of the smallest rectangle containing 
the group, oriented parallel to the direction of the movement. Here we study 
how e depends on the angles (or, aa), and on n. Results are summarized in 
Figure m which shows the average and the extreme values of e over 100 runs as a 
function of the angles {ar, aa), for n = 1, 7, A^ — 1. It is clear that reducing angles 
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from (360° 360°) to (40° 180°) affects the elongation of the group. Wide angles 
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Figure 4: Mean elongation of the group as a function of the sensitivity angles, for 
= 30 and different values of n. The angles are (ar, «&) = (360°— fcl6° 360°— 
fc9°), fc = 0, . . . , 20. Data come from 100 runs. Error bars are the ranges of the 
outcomes. 



(roughly, 360° > > 200° and 360° > > 270°) induce an average elongation 
greater than 1, i.e. the group stretches along the transverse direction. Moreover, 
the range of the outcomes is large, meaning that different initial conditions affect 
significantly the evolution of the system. In the majority of runs the system does 
not reach an equilibrium in a reasonable time, and the simulation is stopped after 
a maximum number of iterations. Conversely, small angles (200° > > 40° and 
270° > tta > 180°) lead to small values of e, with small differences among the 
runs. The steady-state configurations are strongly elongated in the direction of 
motion: in the limit case, we obtain a line as in Figure [5] One can also observe 
a dependence on n: larger values of n show a sharper transition to lines than 
with n ~ 1. The obtained lines are not in general stable: this fact relates to the 
phenomenon of string instability described in the engineering literature [33] : in 
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Figure 5: A line formation, obtained with ai- — 40° Ofa = 180° n = 7, ^ = 10, 
N = 30. The group is moving horizontally from left to right. See text for the 
explanation of the irregularities in the head of the line. 

a line of vehicles which are tracking their forerunners positions, perturbations 
propagate down the line in cascade, leading to instabilities. Control-theoretic 
results presented in [37] state that weakening the interaction forces should reduce 
these negative effects. Consistently, we have observed that damping down to 
zero small repulsion forces improves the stability of the lines. 

Finally, we note that in Figure [5] some "border effects" are visible in the 
head of the line (right side): since in this case we chose n = 7, the animals 
in the front can not interact with a sufficient number of group mates and then 
they do not form a single-file line. 

3.3 Vees 

In this paragraph we show that a restricted frontal repulsion range (ar < 180°) 
induces the formation of V-like patterns. V-like formations have been recently 
obtained in the literature j30] , using an ad hoc model motivated by aerodynam- 
ics considerations. In our model, instead, V-like formations arise as one of the 
anisotropy's effects. Following [20], we adopt a broad definition of V-like forma- 
tions, which includes asymmetric formations (J-like, and echelons) as well. To 
understand the role of anisotropy in the emergence of such patterns, we study 
how the configurations depend on the repulsion angle a^, while we keep fixed 
aa = 36O°0 

^To obtain sharper evidence from simulations, it is useful to strengthen repulsion with 
respect to attraction, taking a value of § larger than for lines: in the following simulations we 
set ^ = 13, and wc also set n = 7 and t)max = 10. 
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In Figure ini we plot the distribution of the angles between the nearest neigh- 
bor and the direction of motion, for four values of ar- A few remarks are in 
order. If a-c = 360° animals do not show any angle preference. If ar = 270°, 
animals show a preference of the front/back positions versus side positions. Fi- 
nally, if Qr < 180°, animals clearly prefer to keep a specific angle, namely 
with respect to their nearest neighbor. 




Figure 6: Distribution of the nearest-neighbor angle for ar = 360°, 270°, 120°, 
60°. Data from 100 runs. 

Second, to obtain more quantitative results about the transition from clus- 
ters to Vees, we introduce an Alignment Index (AI), defined as follows. AI(6') is 
the percentage of individuals whose nearest neighbor is positioned (up to a small 
tolerance Eangic) at a given angle with respect to them. The dependence on 
«!■ of this novel index is shown in Figurc[71 The figure plots both AI(ar/2) and 
AI(30°), computed as the average over 100 runs, with Eangie = 3°. If > 180° 
the alignment index is as low as in a random configuration; if instead a^- < 180° 
the high index confirms the preference for an ^-alignment. Qualitative anal- 
ysis of the obtained configurations confirm these results: one observes that for 
a wide range oi a^, a scenario sets up, in which the animals form (several) V- 
likc formations. Examples are given in Figure [51 Finally, it can be noted from 
Figure [9] that the number of considered neighbors n affects the ability to form 
V-like configurations: if n is too small (71 = 1) or too large (n = N — 1), the 
interesting patterns do not show up. 
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360° 270° 180° 120° 90° 60° 30° 

a 



Figure 7: AI(30°) and AI(ar/2) as functions of a^. Average of 100 runs. The 
reported reference value 6% is the expected value of AI from a random uniform 
distribution. 



^0 -30 -20 -10 
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-19 -18 -17 -16 -15 -14 -13 -12 -11 -10 -9 



Figure 8: V-like formations obtained with Ur = 60° ~ 30. The plot on the 
right is a close-up on one of the Vees. The group is moving horizontally from 
left to right. 



4 Analytical results 

In this section we develop a framework for the analysis of the presented model 
which helps the interpretation of the numerical results. Several analytical tools 
have been developed and used in literature for the analysis of flocking algo- 
rithms. The typical method for their analysis consists in defining a suitable po- 
tential function, called "virtual" or "artificial" potential, whose gradient gives 
the dynamics. Hence a well developed theory on potential systems can be used 
to make a full mathematical analysis, see e.g. [ISJIH]. We have anticipated in 
the introduction that our model has two main features: state- dependent switch- 
ing, and asymmetry of interactions. Let us illustrate them. It is clear from 
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Figure 9: Configurations obtained with n = 1,3,7, N — 1, respectively, and 
ttr — 60° N = 60. The group is moving horizontally from left to right. 

their definition that the attraction and repulsion neighborhoods A" and TZ" 
depend in a discontinuous way on the configuration (i.e. on the positions of 
the animals). As a consequence, the system's evolution switches among a finite 
collection of equations. From a system-theoretical point of view, equation ([3]) 
is a switching system |26j with state-dependent switches. These issues have al- 
ready been taken into account in flocking studies, as many papers consider the 
case of animals which interact when they are closer than a certain threshold. 
See, among others, [38]. The virtual potentials approach can be extended to 
these problems, provided the potential is allowed to be non-differentiable at the 
switching points. However, our model has the distinctive feature that, because 
of the limitation in the number of neighbors and in the shape of the sensitivity 
zones, interactions need not be symmetric. Example giving, for two animals 
the inclusion j e Af does not imply that i G A^. This fact prevents us 
from using a virtual potentials approach, as the operation of differentiating the 
potential does not keep any directionality information. Actually, in literature 
there is no general approach available for systems with directed interactions. 
A partial asymmetry has been taken into account in other works |34| . but the 
given treatment is far for being sufficient for our purposes. In what follows we 
lay down the basics of a theory that we believe is able to catch the specific 
features of our model. 
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Let us start by defining solutions of a switching system in an useful sense. 
Indeed, a differential equation of the form x = f{x) with a discontinuous right- 
hand-side can not have a solution in the classical sense, i.e. a solution which 
is differentiablc. Let us consider the right-hand side of equation ([3]). That 
expression is not defined on the set of degenerate configurations, that is 

S^{xe (R2)^ I 3i,j e {l,...,iV} s.t. X, = Xj}. 

Moreover, it is not always well defined on the set of switching configurations, in 
which two or more agents arc equidistant from another, 

D = {xe (R2)^| 3i,j,fcG{l,...,iV} s.t. \\x,^Xj\\ = \\x,~Xk\\], 

because of the ambiguity in the definition of the "n closest neighbors" . 

To have a proper definition, we shall consider the following differential equa- 
tion 

±{t) = f{x{t)), (4) 

where the flow / : (R^)^ \{DUS) (K^)^ is defined componentwise as in 
equation ([3]), 

Note that / can not be extended with continuity to the set DU S. Hence, 
a solution involving, for instance, two animals equidistant from a third one, 
can not be defined in the classical sense. In what follows, we shall extend 
the definition of the solutions of equation Q to include the set D. For such 
extension, we shall follow the approach in |llj , which requires defining a suitable 
differential inclusion, derived from 

To this goal, let B{y, 6) denote the Euclidean ball of radius 6, centered in y, 
and set 

ICfiy)=r\ n ^{fiBiy,S)\A)}, 

<5>0p(A)=0 

where the operator co denotes closed convex hull, and /i denotes Lebesgue mea- 
sure. The map x : R>o -> (R^)^ is said to be a Filippov solution of the 
system ^ if it is absolutely continuous and it satisfies the differential inclusion 

x{t) G ICf{x{t)) 

for almost any t > 0. 

From now on, we restrict ourselves to the case in which ~ 360°, a,- ~ 360°, 
and n = 1. Hence, each animal interacts just with its closest mate, and A} = TZ\ 
for every i. We make this assumption because, while the analysis for this case 
is simpler than for the general one, still the significant features of switching and 
asymmetry are apparent. Indeed, provided N > 2, the relation of "being the 
closest to" needs not to be symmetrical. Moreover, simulations show for this 
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case the formation of regular structures, as in Figure[3J which can be of intrinsic 
interest. 

In the case we are considering, it is useful to define the set of the closest 
neighbors of a given animal i G {1, . . . , N} as 

closesti(a;) — argmin{||a;i — 

Notice that closesti(a;) may be multivalued when x G D, and let |closesti(a;)| 
denote its cardinality, that is the number of closest neighbors of animal i. 

Given the above definitions, we are able to prove the existence of Filippov 
solutions. 

Theorem 4.1 (Existence) Let n ^ 1. Then, for any initial condition G 
(M^)^ \ S, equation ^ has at least one Filippov solution x, such that x{0) = x^ 
and x{t) e (R2)JV y g ^y^^y ^ > g. 

Proof: We remark that / is piecewise continuous in the following sense. Let 

V = {« e {1, . . . , iV}^ I V* e {1, . . . , N}, V, + i}. (5) 

For any u G V, let us also define the open, possibly empty, set 

£;„ = {.Te (]R2)^\5| ViG {1,...,7V}, < \W ~ x,\\,k i {i^v.^]} 

which is the set of the configurations such that Vi is i's closest neighbor. Note 
that E'l, n £'„ = if u ^ w and the measure of the boundary of each Ey is 
zero. Moreover, denoting by the closure of E^ in the induced topology of 
(R2)W \ it holds that Ut,ev = '^)^ \ In the interior of each "piece" 
Ey, that is if x G £'„, we have that f{x) ~ f^{x), where the function /" is 
defined componentwise as 

f-{j:) = {xy^-x,)-e 



Moreover, each function /" is continuous in Ey. These facts imply that for every 
X G (M^)^ \ 5, the set /C/(x) is bounded, nonempty, closed and convex, and the 
map X ICf{x) is upper semicontinuous. Before we conclude, we need to show 
that, provided x(0) G (R^)^ \ 5", the solution can not reach S either in finite 
time or asymptotically. By contradiction, let 

lim min \\xh{t) ~ Xk{t)\\ = 0, (6) 
t->Lh,ke{i,....N} 

for L G (0, c«]. Note that for each t G [0,i), there exist a pair {h*,k*) G 
{l,...,iV}^, possibly depending on time, such that it attains the minimum 
in ([SJ, and thus h* ~ closest^* (.T(t)) and k* = closest^. (.T(t)). By continuity of 
the solution, it exists to G (0, L) and e G (0, ^) such that for t G [to, L), we have 
I \xh-' (t) — Xk* {t)\ \ < £ < This implies that repulsion is larger than attraction. 
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and then the animals are moving away from each other, that is \ \xh* (t) — Xk* (t) \ \ 
is increasing in t. This contradicts equation ([6]). 

At this point, we can apply [TTl §7 Theorem 1], and obtain that the differ- 
ential inclusion i{t) € lCf{x(t)) has at least one solution x{t), for alH > and 
for any initial condition a:° G (M^)^ \S. □ 

Remark 4.2 The above proof can be extended to the case n > 1 and < 360° 
ar < 360° modulo a suitable redefinition of the "pieces" E^, in order to account 
for the more complex neighborhood relationships: the notational setup would 
be cumbersome, and wc do not detail it. 

Loosely speaking, we would expect that a configuration having the closest 
neighbor(s) at distance ^ for all animals, as in the lattice configuration of Fig- 
ure[3a]would be an equilibrium configuration. Indeed, each animal i is driven by 
the attraction-repulsion force component fi towards keeping a distance ^ from 
its neighbors. The following result technically clarifies this intuition. Before the 
statement, we need to define a configuration x* to be a Filippov equilibrium of 
/ when e K f{x*). 

Proposition 4.3 (Equilibria) Let x* e {M.'^)'^ \ S. If for all i e {1, ... ,7V} 
and for all k e closest; (a;*), it holds \\x* — x'^\\ = ^, then x* is a Filippov 
equilibrium for the system and moreover 1 < |closcsti(a:;*)| < 6, for all 
i€{l,...,N}. 

Proof: Wc consider closest; (x*) for every i € {!,..., A^}, and distin- 
guish the cases in which their cardinalities are equal to, or larger than 1. If 
|closesti(a;*)| = 1 for all i e {1,...,A^}, then / is smooth at x* , and hence 
ICf{x*) = {fix*)} = {0}. If instead there exists h G {l,...,iV} such that 
I closest?! (x*) I > 1, then let us consider the set V defined in ([5]), and take v gV 
such that Vi G closest; (.t*), for every i G {!,..., N}. Since x* is an accumula- 
tion point of El,, let us consider a sequence {cc'jigN C E^, such that x^ — > x* as 
I — >■ -l-oo. Hence, /(a;') — > as Z — > oo, and this implies, by the definition of the 
differential inclusion, that G ICf{x*). 

By definition, | closest; (a;*)| > 1. The fact that |closest;(a;*)| < 6 can be 
shown by contradiction. Let |closesti(a;*)| = m, with m > 6, for some i. Then 
there are m points of M^, representing positions, which belong to a circle of 
radius ^ centered in x* . But then the distance between two of them has to be 
less than ^, which is a contradiction. □ 

Hexagonal lattice configurations, in which animals are at distance ^ from 
their closest ncighbor(s), arc observed in simulations in the case aa = = 360° 
and n = 1, as reported in Figure [5al Indeed, Proposition 14. 31 shows that such a 
lattice is an equilibrium: loosely speaking, we can say that it is the most closely 
packed configuration among the equilibria pointed out by this result. Notice 
moreover that such lattice equilibria are actually switching configurations, be- 
longing to the set D: this gives an a posteriori justification of the effort that we 
have done for a careful extension of the solutions to this set. 
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5 Discussion 



Our agent-based model has been conceived from assumptions which are widely 
accepted from the biological point of view, and it shows the onset of group 
structures and patterns which are observed in nature. In some sense, our model 
can be seen as a generalization of other similar models, since in the special 
case (ttr, tta) = (360° 360°) we recover results which are by now consolidated in 
literature. The novelty resides in that by modifying a small set of key param- 
eters, we obtain other patterns, which are experimentally observed in animal 
groups. Our results suggest that apparently large differences in group patterns 
may arise just from differences in the attraction and repulsion sensitivity zones. 
In this perspective, we would like to point out that the group structure is not 
only a function of the species, but also of the external conditions. For instance, 
surf scoters |27| and other animals can form either clusters or lines, depending 
on the environmental conditions. This suggests that and are not only 
animal-dependent, but they can also vary depending on the type of motion, en- 
vironmental conditions, presence of predators, and aim of the displacement. In 
this paper, the choices of the angles were made to explore the model's properties 
and they do not come from experimental data. Nevertheless, the correspondence 
of the simulated patters with those experimentally observed suggests that fu- 
ture research should distinguish between the (variable) sensitivity zones and the 
(fixed) visual field, as in [HI [30] . 

Now we review the results obtained in the previous sections. The goal is 
to propose some biological insights, especially about the role of anisotropy in 
the interactions. As we did presenting the simulations, we distinguish three 
cases: isotropic attraction and repulsion; anisotropic attraction and repulsion; 
isotropic attraction and anisotropic repulsion. 

In the first case, with completely isotropic interactions, our model produces 
clusters of individuals. Clusters are common for small birds (e.g., starlings [3]) 
and for fish, whose sensing abilities allow an almost complete perception around 
them. Moreover, these animals often change the direction of motion following 
complex trajectories, and this makes useful to keep under control all the space 
around them. By means of isotropic interactions, they are able to change rapidly 
the direction of motion of the whole group, with no need to significantly rear- 
range the shape of the group. This would not be possible in less symmetric 
formations like lines or Vees. 

Instead, when interactions are not isotropic, also the obtained formations are 
not isotropic, and we observe that their onset depends on spontaneous leader- 
following mechanisms. When both repulsion and attraction are focused in front, 
the leader- following mechanism induces the formation of lines. Lines are com- 
monly observed in slow-moving animals as lobsters, elephants and penguins. 
The small repulsion angle ai- appears to be crucial to obtain such a pattern (see 
Figure [4]). This could be related to the fact that such animals keep a steady 
direction of motion, and, once the line is established, repulsion needs to be ac- 
tive only against the forerunners. The choice — 180° can also be significant 
from the biological point of view. Indeed, this would mean that individuals pay 
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more attention to the group mates in front, while they do not respond to what 
happens behind them. For example, they would not perceive a disconnection of 
the group. Hence, a model assuming = 180° seems suitable either for animals 
with a restricted frontal visual field or for animals which are not particularly 
interested in the cohesion of the whole group. 

A different leader-following mechanism induces V-likc formations when at- 
traction is isotropic and repulsion is restricted to the group mates in front. Our 
simulations reproduce several significant features of natural V-like formations, 
described in the experimental literature [T3J [501 HZ] . First, echelons and J-like 
formations are as common as perfect Vees. In our model all these formations 
appear, and actually the behavior of single individuals does not depend on the 
global shape of the formation. Second, V-like formations are not stable, but 
rather they often disband and quickly reform. Third, in our model each Vee is 
made of a limited number of individuals (see Figure [8]), independently of the 
total number of animals N. In other words, increasing N leads to the formation 
of a larger number of V-like groups, but not to larger groups. This fact has 
been experimentally observed [17] and explained [32| resorting to the argument 
of string instability, that we introduced in Section 13.21 Our results should be 
compared to those in [30], where authors propose an ad-hoc formation algo- 
rithm based on aerodynamic arguments in which the number of V-like groups 
is constant for increasing A^. 

The function of V-like formation has been greatly investigated (see for in- 
stance [17l [To] [35] [421 [32] and references therein). Two hypothesis are the most 
considered: aerodynamic advantage, and visual communication advantage. The 
former is based on the fact that each flying individual creates an upwash region 
behind it, just off the tips of its wings, so that another individual can benefit 
placing itself in that region. The latter is instead based on the fact that flying 
in a skewed position with respect to the bird in front is useful to avoid collisions, 
and allows an unhindered visual communication with all the group mates |10| . 
Our results appear to support the hypothesis of visual communication advan- 
tage, since V-like formations are not obtained imposing individuals to stay in 
upwash regions. Instead, they are obtained just from frontal repulsion and, for 
the sake of group cohesion, isotropic attraction. 

Motivated by some literature [Tl [lH [2] , we have also investigated the depen- 
dence of the configurations on n, the number of interacting neighbors. We have 
seen that cluster configurations can be obtained with any n, but the value of 
n influences the internal structure of the cluster. Moreover, it appears that an 
intermediate number of neighbors is preferable to obtain lines or V-like forma- 
tions. These results seem to confirm the reasonableness of topological interac- 
tions with a limited number of neighbors. However, this point deserves further 
investigations. 

Future research 

We are keen on developing our research in three respects. First, it should be 
noted that although the presented model is bidimensional, it can readily be 



18 



extended to a three-dimensional one. We are planning to study a 3D version of 
our model since the novel set of simulations might show interesting phenomena. 
We are especially interested in studying the effect of the anisotropy in 3D V- 
like formations, because the skewed formations described in |20j are inherently 
t hree- dimensional . 

Second, we want to further develop the theoretical analysis of the model, 
in order to include a thorough description of the equilibria of ([2]), and their 
stability analysis in the suitable switching systems framework. Moreover, we 
would be interested in a variational interpretation of the model. Indeed, while 
it is apparent that each animal is trying to minimize a "private" potential de- 
pending on its neighbors, it is unclear whether and when this would result in a 
configuration minimizing some global objective function. 

Third, we believe that the simplicity of the interaction rules we have pro- 
posed, and the limited number of group mates to be kept into account, can be 
interesting features for engineering applications. In particular, we will investi- 
gate the design of robust algorithms for environmental deployment of robots or 
sensors and the formation cruising of unmanned vehicles. 
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